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Abstract: We study a fixed point formalisation of the well known analysis of Bianchi. We 
provide a significant simplification and generalisation of the analysis. In this more general 
framework, the fixed point solution and performance measures resulting from it are studied. 
Uniqueness of the fixed point is established. Simple and general throughput formulas are 
provided. It is shown that the throughput of any flow will be bounded by the one with the 
smallest transmission rate. The aggregate throughput is bounded by the reciprocal of the 
harmonic mean of the transmission rates. In an asymptotic regime with a large number 
of nodes, explicit formulas for the collision probability, the aggregate attempt rate and the 
aggregate throughput are provided. The results from the analysis are compared with ns2 
simulations, and also with an exact Markov model of the back-off process. It is shown how 
the saturated network analysis can be used to obtain TCP transfer throughputs in some 
cases. 
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Nouvelles perspectives d'une analyse de point fixe d'une 
seule cellule de WLANs IEEE 802.11 


Résumé : Nous étudions une formulation d'un point fixe de l'analyse connue de Bianchi. 
Nous proposons une simplification considérable et généralisons l'analyse. Nous étudions 
dans ce cadre plus général, la solution du probléme de point fixe ainsi que les mesures de 
performances qui en découlent. Nous démontrons l'unicité de la solution du probléme de 
point fixe. Des formules simples et générales du débit sont obtenues. Nous montrons que 
le débit d'un flux quelconque est borné par le taux de transmission le plus petit. Le débit 
aggrégé est borné par l'inverse de la moyenne harmonique des taux de transmissions. Dans 
le régime asymptotique d'un grand nombre de noeuds (mobiles), nous obtenons des formules 
explicites pour la probabilité de collision, le taux aggregé de tentatives de transmission et 
du débit aggrégé. Les résultats de l'analyse sont comparés aux simulations par ns2 et aussi 
par une approche markovienne exacte du processus de back-off. Nous montrons comment 
l'analyse du réseau saturé peut étre utilisé pour fournir le débit de connexions TCP dans 
certains cas. 


Mots-clés : IEEES02.11, Évaluation de performances, Equation point fixe, CSMA/CA, 
Réseaux sans fil 
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1 Introduction 


We consider wireless local area networks that use the IEEE 802.11 standards for physical 
transmission and medium access control (MAC) (see [8] for a concise description of these 
standards and some planned extensions). In particular, we are concerned in this paper 
with the network architectures depicted in Figure 1. There are several nodes within such 
a distance of each other that only one transmission can be sustained at any point of time. 
We call these single-cell networks (see Figure 1). Our discussion covers networks of both 
the types shown in the figure, i.e., ad hoc networks in which there is no central access point 
(AP) through which all traffic must pass, or infrastructure networks in which an AP acts 
as a conduit between the wireless network and a wired "infrastructure." In addition, in this 
paper our analysis is limited to the situation in which all nodes use the RTS/CTS based 
distributed coordination function (DCF) without the QoS extensions (as in IEEE 802.11e) 
(but see [9] for our extensions of the work in the present paper). 

There are several nodes, each of which may have several physical connections or asso- 
ciations with several other nodes. On each such connection the sustainable physical trans- 
mission rate may be different. The physical layer definitions in the IEEE 802.11 standard 
support multiple transmission rates, and the standard includes mechanisms that permit 
nodes to agree on a transmission rate that performs well between them. Between each such 
pair of nodes there are flows whose throughput performance we are concerned with. It is 
assumed throughout this paper that all flows are infinitely back-logged at their transmitters; 
i.e., there are always packets to transmit when a node gets a chance to do so. 

In such a scenario, we are interested in obtaining quantitative formulas and qualitative 
insights via a stochastic analysis of the way that the IEEE 802.11 CSMA/CA protocol 


Figure 1: Single cell IEEE 802.11 networks; the one on the left is in the ad hoc mode, and 
the one on the right is in the infrastructure mode, with the AP serving as the access to a 
wired network. In the left panel C; are the achievable physical transmission rates on the 
corresponding connections. In the right panel Cu; refer to “up-link” transmission rates, and 
Ca, refer to “down-link” transmission rates. In the IEEE 802.11 standard a receiver and a 
transmitter agree on a common physical transmission rate. 
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allocates the wireless medium to the node transmitters. Our approach is to begin with a 
key approximation made by Bianchi [3]. This leads to a fixed point equation, which can 
be expected to characterise the operating points of the system. This fixed point equation 
is our point of departure. We simplify and generalise the analysis leading to the fixed 
point equation. We then establish a simple, and practically appealing, condition for the 
uniqueness of the fixed point in this more general framework. Some simple observations 
lead to throughput formulas for the overall network and for the individual flows. These 
formulas allow us to recover the well known observation that the slowest transmission rate 
dominates the throughput performance. We also analyse the fixed point in the asymptotic 
regime of a large number of nodes and find explicit formulas for the collision probability, 
the channel access rate and the network throughput. A key parameter in the protocol is 
the back-off multiplier, whose default value in the IEEE 802.11 MAC standard is 2; we 
numerically study how the back-off multiplier affects the throughput. 

We provide ns2 simulation results for the collision probabilities and throughputs and 
compare these with results obtained from the fixed point analysis. We also provide results 
from an exact Markov chain model for the back-off process and also compare these results 
with those from the fixed point analysis. 

As already pointed out, the above described modeling assumes that there are always 
packets backlogged on every connection. Such a saturation assumption is a common simpli- 
fication and is useful in the following ways. In some situations it has been formally proved 
(see, for example, [6]) that the saturation throughput provides a sufficient condition for 
stability of the queues; ie., if at each queue the arrival rate is less than the saturation 
throughput then the queues will have a proper, joint stationary distribution. In this paper 
we also apply the saturation throughput analysis to provide an analysis for TCP controlled 
file transfer throughputs in certain local area network scenarios. 

Modeling and analysis of wireless networks is extremely important as, in addition to 
providing performance results, it helps to understand the strengths and limitations of the 
evolving protocols and standards. The role of such modeling in interpreting experimental 
results and in checking simulation outputs cannot be over-emphasised. The most popular 
model for IEEE 802.11 networks, and one that has led to many applications and extensions, 
is the one reported in [3]. Another analysis, that also incorporates the feature of adapting 
the back-off parameters, has been reported in [4]. The recent paper [2] is one of the many 
that have reported a throughput “anomaly” in IEEE 802.11 networks; i.e., if the network 
has low speed connections, even the high speed connections experience throughput no better 
than what is obtained by the low speed connections. 

The paper is organised as follows. In Section 2 we provide the key observation and 
approximation on which the analysis is based. In Section 3 we analyse the back-off process 
in à fairly general setting. The fixed-point equation is provided in Section 4 and analysed 
in Section 5. In Section 6 the throughput formulas are provided. The asymptotic analysis 
is developed in Section 7. The longer proofs are in the Appendix. 
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Figure 2: The evolution of the back-off periods and channel activity for four nodes. It can be 
seen that back-offs are interrupted by channel activity, i.e., packet transmissions and RTS 
collisions. 


2 A Key Observation and an Approximation 


2.1 Sufficient to Analyse the Back-Off Process 


We begin by extracting from a description of the system the key modeling abstractions that 
will allow us to develop the analysis. Figure 2 shows the evolution of the system for 4 nodes; 
shown are the back-offs, the transmissions and collisions. In the IEEE 802.11 standard, the 
back-off durations are in multiples of a standardised time interval called a slot (e.g., 20 us in 
IEEE 802.11b). However, this discrete nature of the back-offs does not affect the following 
argument. When a node completes its back-off (for example, node 1 is the first to complete 
its back-off in Figure 2), it seeks a reservation of the channel by sending an RTS packet. If no 
other node completes its back-off before hearing this transmission then the RTS effectively 
reserves the channel for the first node. There follows a CTS from the intended recipient of 
the RTS, and then there follows a packet transmission and a MAC level ACK. This ends 
the reservation period and the node that transmitted the packet samples a new back-off 
interval. Note that we assume throughout that nodes always have packets to transmit; i.e., 
all the transmission queues are saturated. 

If the RTS collides with that of another node, then after fully transmitting their RTSs 
each node waits for a time interval DIFS before returning to the back-off state. For example, 
in Figure 2 nodes 2 and 3 collide after the first two attempts (by nodes 1 and 3, respectively) 
are successful. The other nodes, not involved in the collision, listen to the channel activity 
until the end of the RTS transmissions, and then also start their DIFS timers. Thus after 
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Figure 3: After removing the channel activity from Figure 2 only the back-offs remain. 
These determine the scheduling of channel access. At the bottom is shown the aggregate 
attempt process on the channel, with three successes and one collision. 


a collision all nodes resume their back-off phases after an amount of time equal to the 
transmission time of an RTS plus a DIFS. 

If attempts to send the packet at the head-of-the-line (HOL) meet with several successive 
failures, this packet is discarded. By our assumption of saturated queues, there is always 
another packet waiting to be sent by the upper layers: either the same packet or the next 
one in line. 

We see from the figure that when any node has reserved the channel or whenever there 
is a collision, all other nodes freeze their back-off timers. We also notice that the evolution 
of the channel activity after an attempt is deterministic. It is either the time taken for a 
transmission or for a collision. If there is a transmission then the time depends on which 
node captures the channel. The latter dependence comes about because the transmission 
time of a packet depends on the transmission rate and hence on the transmitting node. 

Since all nodes freeze their back-offs during channel activity, the total time spent in 
back-off up to any time t, is the same for every node. With this observation, let us now look 
at Figure 3 which shows the back-offs of Figure 2 with the channel activity removed. Thus 
in this picture “time” is just the cumulative back-off time at each node. In the IEEE 802.11 
standard the back-offs are multiples of the slot time. A success occurs if a single back-off 
ends at a slot boundary, and a collision occurs when two or more back-offs end at a slot 
boundary. The nodes could have different back-off parameters (the mean back-off intervals, 
how these are varied in response to collisions and successes, and the number of retries of 
a packet). It is clear, however, that the (random) sequence in which the nodes seek turns 
to access the channel and whether or not each such attempt succeeds depends only on the 
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back-off process shown in Figure 3. It is therefore sufficient to analyse the back-off process 
in order to understand the channel allocation process. The saturation assumption is crucial 
here since, with this assumption, we do not have to take care of any external packet arrivals 
that may occur during channel activity periods. 

Thus, in summary, we can delete the channel activity periods, and we are left with a 
"conditional time" which we will call back-off time. We will analyse the back-off process 
conditioned on being in back-off time. It will then be shown how this analysis can be used 
to yield the desired performance measures over all time. 


2.2 A Key Approximation 


Throughout the rest of the paper we assume that all the nodes use the same back-off pa- 
rameters. Hence the back-off process shown in Figure 3 is symmetric over the nodes. We 
call this the homogeneous case to distinguish it from the nonhomogeneous case in 
which different nodes may use different back-off parameters, as, for example, proposed in 
the IEEE 802.11e standard (see [8]). 

In Figure 3 we also show the aggregate sequence of successes and collisions. In general, 

this is a complex process, and it is also clear that the success and collision processes of the 
various nodes are coupled and strongly correlated. In Section 8 we will describe an exact 
Markov chain model for the joint backoff process of the nodes, but this model is analytically 
intractable. The following key approximation is made in [3]. 
The Decoupling Approximation: Let Ó denote the long run average back-off rate (in 
back-off time) for each node. By the fact that all nodes use the same back-off parameters, 
and by symmetry, it is assumed that all nodes achieve the same value of 8. Let there be n 
contending transmitters, and consider a given node. The key approximation is to assume 
that the aggregate attempt process of the other (n — 1) nodes is independent of the back-off 
process of the given node. In IEEE 802.11 the back-off evolves over slots, hence a discrete 
time model (embedded at slot boundaries) can be adopted. Then the approximation says 
that if the attempt rate per node is 3 attempts per slot (0 € 8 < 1), then from the point 
of view of the given node the number of attempts by the other nodes in successive slots 
are i.i.d. Binomial random variables with parameters (n — 1) and B. It might be expected 
that such a decoupling approximation should work well when there is a large number of 
transmitters accessing the channel. 


3 Analysis of the Back-Off Process 


We generalise the back-off behaviour of the nodes, and define the following back-off param- 
eters. 


K :— At the (K + 1)th attempt either the packet succeeds or is discarded 


b, :— The mean back-off duration (in slots) at the kth attempt for a packet, 0 € k < K 
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Figure 4: Evolution of the back-offs of a node. Each attempted packet starts a new back-off 
“cycle”? 


Since we are limiting ourselves to the homogeneous case, these parameters are the same for 
all the nodes. 

In Figure 4 we show the evolution of the back-off process for a single node. There are 
R; attempts until success for the jth packet (no case of a discarded packet is shown in this 
diagram), and the sequence of back-offs for the jth packet is B. O<i<R;-—1. Thus the 
total back-off for the jth packet is given by 


with e(a) = b;. We observe that the sequence X;, j > 1, are renewal life times. Hence, 


viewing the number of attempts R; for the jth packet as a "reward" associated with the 
renewal cycle of length X;, we obtain from the renewal reward theorem that the back-off 
rate is given by as. 


Now let y be the collision probability seen by a node, i.e., 
y := Pr (an attempt by a node fails because of a collision) 


Since the back-off behaviour of all the nodes is the same, the collision probability is the same 
for all the nodes. By the approximation made in Section 2, the successive collision events 
are independent. It is then easily seen that 

E(R) = 14+7+7---+7*% 

E(X) bo + ybi + 253 + 4 y b + yÉ b 


which yields the following formula for the attempt rate for a given collision probability 7 


l+y+7---+7* 


) bo + ybi + 3303 +o be + ROK 
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Figure 5: Evolution of the back-off stage (Zi) and the residual back-off time (Y (t)) for the 
case in which the back-offs are continuous variables. 


(1) 
Note that, since the back-off times are in slots, the attempt rate G(y) is in attempts per slot. 
Remarks 3.1 


1. Note that the distribution of the back-off durations does not matter. 


2. Also, observe that the above analysis remains unchanged whether the back-off distri- 
butions are discrete (i.e., the back-offs evolve over slots) or are continuous. 

3. In the back-off model considered in [3] K = oo; further, there is an m > 1 such that 
b, = (Zeupaai) slots, for 0 € k € m — 1, and by = (2 Clans) slots, for k > m. 


Here C Wy, is a positive integer (2° in the IEEE 802.11 standard). Substituting these 
into the expression for G(y) in Equation 1 yields 


2(1— 29) 


G0) = Cas € D + 30Was d = 037] 


attempts per slot, which is the same as in the paper. Note that the + alternatives 
arise depending on whether we take the back-off to be uniformly distributed over 
[1,2,---, CW] or over [0, 1, ---, CW — 1]. Evidently, the uniform distribution of back- 
off durations plays no role in the final results in [3]. 


4. À more detailed evolution of the back-off process in Figure 4 is shown in Figure 5, 
where at each time ¢ the residual back-off duration Y(t) is also shown. The process Zp 
is the back-off stage the node is in. Thus if K = 7, Zi = 3, and Y(t) = 5, then after 5 
time units the current back-off ends. If there is a collision, Z, = 4 and a back-off with 
mean b4 is sampled from the specified back-off distribution (uniform in the standard). 
If Zk = 7 then at the end of the current back-off, irrespective of whether there is a 
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collision or a success, the next back-off has mean bo, and is sampled from the specified 
distribution. It is clear that the process (Zx, Y (t)) is Markov regenerative, with Zi 
being the embedded Markov process. It is, however, well known that it is unnecessary 
to analyse (Zi, Y (t)) to obtain event rates. The process Z; and the successive back-off 
intervals constitute a Markov renewal process. It is also well known that for Markov 
renewal processes event rates and time probabilities are insensitive to distributions of 
life-times. It should thus be clear why one can directly obtain the formulas above 
without needing to go through the analysis of the Markov chain in [3], and also why 
the results are insensitive to the back-off distribution. 


4 The Fixed Point Equation 


Focusing on the back-off and attempt process of a node, and being given the collision prob- 
ability y the attempt rate is provided by G(y) in Equation 1. It is important to recall that 
in the present discussion all rates are conditioned on being in the back-off periods. Later 
we will see how to incorporate the channel activity periods. Now if all nodes have the same 
back-off parameters, they will all see the same average collision probability, y, and hence 
will have the same attempt rate. If the attempt rate (or probability) of each node per slot is 
8,0 € 8 € 1, then, conditioning on an attempt of the given node, the probability of this at- 
tempt experiencing a collision is the probability that any of the other nodes attempts in the 
same slot. Under the decoupling approximation, the number of attempts made by the other 
nodes is binomially distributed with parameters @ and n — 1. Under the approximation, the 
number of attempts in successive form an i.i.d. sequence. The probability of collision of an 
attempt by a node is given by 


DB) := DES (2) 


We will show later in the paper that under a certain asymptotic regime the aggregate attempt 
rate nj converges to a positive value as n — oo. Then (motivated by the binomial to Poisson 
convergence theorem) for a large number of nodes, it is reasonable to model the attempt 
process of the other nodes (with respect to a given node) as a sequence of i.i.d. batches (at 
slot boundaries) with the batch distribution being Poisson with mean (n — 1). The collision 
probability under this model is then clearly given by 


T(8) :2 1-e-(-U8 (3) 


It is now natural to expect that the equilibrium behaviour of the system will be charac- 
terised by the solutions of the following fixed point equation 


y = Y(G(y) (4) 


If this equation can be solved it will yield the collision probability, from which the attempt 
rate can be obtained using Equation 1. We will see in Section 6 that throughputs can be 
obtained once these quantities are determined. 

When faced with a fixed point formulation, several questions need to be asked. 
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1. Existence and uniqueness: While fixed points can usually be shown to exist, they 


may not be unique. Existence of multiple fixed points may suggest multistability, i.e., 
the system may exist in one of several operating states for long periods of time. Often 
at least one of these operating states has poor performance; see, for example, [1] and 
[5] for an example of bistability in circuit switched networks with adaptive routing. 
We will provide conditions for the uniqueness of the fixed point. 


. Computation of the fixed point: If the function l'(G(*)) is a contraction then the 


fixed point iteration, i.e., 
^1 = T(G(o)) 


will converge starting from any yo € [0,1]. As is often the case, we do not have the 
contractive property. The following relaxed fixed point iteration may work with an 
appropriate choice of o € (0,1) 


Ve = (1 — o)T(G(ys)) + or 


The choice of a that leads to convergence will, in general, depend on parameters of 
the problem. We will examine this technique for the special case in which K — oo in 
Section 7.1. For finite K, the numerical results provided in the paper are obtained by 
using the solve() function in MATLAB. 


. Convergence of the system to the fixed point in some sense: Even when a 


fixed point exists and is unique it is not immediately clear in what way the steady state 
system performance will be related to the quantities obtained from the fixed point. 
One needs to show that the system dynamics converge in some sense and appropriate 
quantities in the limit satisfy the fixed point equation. Such an analysis may also lead 
to a justification of the assumptions that lead to the fixed point equation (e.g., the 
Poisson and independence assumption in the present paper). See [5] for an example 
of such an analysis. À related analysis in the context of packet switching systems has 
been done in [10]. We will defer this issue for our future investigations. In the present 
paper we proceed assuming that the steady state system behaviour will correspond to 
the unique fixed point. 


. Analysis of the fixed point to understand system performance: The fixed 


point equation can be analysed to characterise the fixed point and thus obtain quanti- 
tative and qualitative results about steady state system performance. The major part 
of this paper is devoted to such an analysis. 


Analysis of the Fixed Point Problem 


Since l'(G(y4)) is a composition of continuous functions it is continuous. We thus have a 
continuous mapping from [0, 1] to [0,1]. Hence by Brouwer's fixed point theorem there exists 
a fixed point in [0, 1]. We next turn to uniqueness. 
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Lemma 5.1 G(»y) is non-increasing in y if by, k > 0, is a non-decreasing sequence. 
Proof: See the Appendix. 


Theorem 5.1 T'(G(*)) : [0,1] — [0,1], has a unique fired point if bk, k > 0, is a nonde- 
creasing sequence. 


Proof: Since T'(B) is non-decreasing in 8 and, by Lemma 5.1, G(y) is non-increasing in 
q, it follows that [(G(y)) is non-increasing in y. The fixed point must therefore be unique, 
since multiple fixed points will lead to a contradiction to the non-increasing property of 
T(G(5). m 
Remark: We observe that in the IEEE 802.11 standard the sequence b; is non-decreasing. 
Hence for the practical system there will be a unique fixed point. 


5.1 Examples of the Fixed Point Solution 


In Figure 6, we show plots of l'(G(»)) vs. y for several parameters. The back-off is expanded 
multiplicatively with multiplier p — 2, as in the IEEE 802.11 standard. In the plot on the 
top we use the value K = 7 which corresponds to the retry limit in the standard. In both the 
plots the initial mean back-off b is 16 slots. The intersection of these plots with the *y—x" 
line corresponds to the fixed point. We see that the collision probability increases with an 
increasing number of nodes. For n > 30, with K = 7, the collision probability is larger than 
with K = 100. The collision probability for n < 20 does not appear to be sensitive to K for 
K > 7. 

As an illustration, in Figure 7 we show fixed points for some “bad” back-off adaptations. 
The retry limit is 2 in both the plots, the number of nodes is 20. In the plot on the top 
a collision causes the back-off to be reduced by a factor of .01. We see that there is a 
fixed point very near to 1, i.e., a collision probability of 1, implying very poor performance. 
In the plot on the bottom the initial back-off is also adapted inversely with the collision 
probability (i.e., bo = zu and in addition the back-off is reduced by a factor of .01 when 
there is a collision. There are now two fixed points, one with a collision probability close to 
zero, implying an idle network, and the other with a collision probability close to 1 implying 
that all attempts are wasted. With this, obviously bad, back-off adaptation one can expect 
bi-stability between two very poorly performing regimes. These contrived examples have 
been included to serve as a caution that attempts to modify the back-off mechanism have 
to be made with care. 


5.2 Comparison with an ns2 Simulation 


It was reported in [3] that the fixed point analysis works well for IEEE 802.11 parameters. 
In Figure 8 we demonstrate this by plotting the collision probability obtained from the fixed 
point method and from an ns2 simulation. 

In all the ns2 simulations presented in this paper we have used ns2 version 2.26. The 
bugs present in the IEEE 802.11 code were patched by using an updated version of the code 
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Figure 6: Plots of [(G(7)) vs. y for two values of K (7 (top) and 100 (bottom)), bo = 16 slots, 
and multiplicatively increasing b, with multiplier p = 2. For each K, plots are shown for 
number of nodes n — 2, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100. 


taken from the ns2 snapshot dated January 5, 2004. Static routing was implemented by 
using NOAH code (dated November 2003), downloaded from the web site of J. Widmer, 
EPFL, (http:/ /icapeople.epfl.ch/widmer/uwb/ns-2 /noah/index.html). 

As can be seen, the fixed point analysis provides a good approximation for a wide range 
of values of the number nodes. 


6 Calculating Throughputs 


Having obtained the collision probability y and the attempt rate per node 8 (during back-off 
periods), we now turn to the problem of calculating the throughputs. In order to do this we 
first make two key observations. 

'The first is demonstrated by Figure 9. Because of the i.i.d. batch binomial assumption on 
the aggregate attempt process, the instants at which a successful transmission or a collision 
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collision probability (y) 
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0 02 04 06 08 1 
collision probability (y) 
Figure 7: Plots of [(G(7)) vs. y for K = 2, and multiplicatively increasing by with multiplier 


p — .01. In the plots on the top bo = 10 slots, and in the plots on the bottom bo = E slots. 
For each case, plots are shown for number of nodes n = 20. 


ends are renewal instants. Each such instant is followed by a time until the next attempt, 
followed by a collision or a success, and so on. 

The second observation is that since all the nodes follow the same back-off process, each 
node has an equal probability of winning the allocation “race.” With this in mind we can 
now discard the back-off times and focus only on the times when an attempt is made and on 
the intervening channel activity. À successful attempt leads to the channel being allocated 
to one of the n contending nodes with equal probability. Hence in a saturated system, in 
order to compute the amount of time the channel will be allocated to a node, we only need 
to know the identity of the packet that will be found at the head-of-the-line if the channel 
is allocated to the node. 

Consider the model shown in Figure 10. The nodes are visited in random order with 
equal probability. Each node receives an open loop stream of packets. There are m; streams 


INRIA 


Fixed point analysis of IEEES02.11 15 


0.6+ 


FPA 
— ns2 
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Figure 8: Plot of collision probability versus number of nodes. Comparison of collision 
probability (y) obtained from an ns2 simulation (plot labeled ns2), and the fixed point 
analysis (plot labeled FP). 95% confidence intervals are shown for the values obtained from 
the ns2 simulation. In the ns2 simulation the default IEEE 802.11 parameters are used: see 
Remarks 3.1, paragraph 3, with K = 6 and m = 5; the data rate is 11 Mbps and the control 
packet rate is 2 Mbps. 


being handled by node i. These are indexed by 1 € j € m,; these would represent m; flows 
from node i to some of the other nodes. We can thus use the term “flow (i, j)”. 

By “open-loop” we mean that packets arrive to the node and have to be delivered; there 
are no acknowledgement and flow control as in TCP controlled traffic. A fraction p;,; of the 


back- back- 
off Success off 


collision 
o 
£5 
O 
T 


renewal instants 


Figure 9: The aggregate process of back-offs and channel activity 
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5 of visit, with 
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p. 
n,j n 


packet "routing" 
probabilities 


Figure 10: The n transmitters are served in random order with equal probability for each 
node. 


packets at node à belong to stream j,1 < j € m;. Since the node is saturated there is always 
a packet at the head-of-the-line when the channel is allocated to any node, and p;,; is the 
probability that the packet is from flow j. Let us define the packet length of flow (i, j) to 
be Lj; and the physical transmission rate for flow (i, j) to be C;,; bits per slot. 

In addition, we define, 


T, := is the fixed overhead with a packet transmission in slots (e.g., IEEE 802.11b: T, = 90.8 
slots) 


T, := is the fixed overhead for an RTS collision in slots (e.g., IEEE 802.11b: T; = 17 slots) 


The above two observations, the traffic model described above, and the parameters listed 
above, lead immediately to the expression in Equation 5 for the saturation throughput of 
flow (i,j) (in bits per slot) given the collision probability y and the per node attempt rate 


0,5(B) = (5) 
8ü-8y- 
ERE D pijL 


The formula follows from the renewal reward theorem. The mean renewal time (see 
Figure 9) is the mean time until an attempt, pus the mean time for channel activity; i.e., 
a transmission or a collision. The term OT is the mean time until an attempt and 
assumes that the aggregate attempt process is binomial. When there is an attempt the 
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channel is allocated to node à (with probability fue). else there is a collision, for 


which the channel will be busy for the time T,. If the channel reservation succeeds, then 


the head-of-the-line packet at node i is of flow (i, j) with probability p; j, and transmitting 
: : Li,j : : 1—8)^-1 

this takes the time ms + To. The mean reward during the cycle is re Pig Dij; thus 

yielding the displayed expression. Canceling the term 1 — (1— 8)", the formula simplifies to 

the expression in Equation 6. 


0,5(8) = (6) 
Bü — p) pig Lis 
1+ Xx (60 - 87 (Oa vin Bet) +To)) + (1 Q = 8)" = npl = BYP) T.) 


6.1 Low Speed Transmitters Bound All Throughputs 


It has been observed (see, for example, [2]) that when there are several flows with different 
physical transmission rates then the throughput of all the flows is bounded by the slowest 
transmission rate. We can examine this observation using Equation 6. 

If 2 nodes 71 and i» are such that for some j1, 1 € jı € Mu, and j2, 1 < jo < M2 


Pir Li, Ji T Pisa Li, „j2 
then it follows from Equation 6 that 
Gi, Y (^; B) < min{Ci E Ci zd 


and 
0i, Ja (7, B) < min(Ci, J) Ci, gat 


i.e., the flow with the lower physical rate will bound the throughput of both. 
Remark: The above analysis points to an important observation. Suppose we are interested 
in achieving flow throughputs that are proportional to their physical link rates; i.e., 0; ; = 
VCi, for some v. It has been suggested in previous literature that this can be achieved 
by appropriately choosing the packet lengths. We notice from Equation 6 that the desired 
throughput proportionality can be achieved only be making L;,; proportional to sai , which 
requires knowledge of the p; js, which may not be practicable. | 

Let us now consider a simpler situation with n nodes each being the transmitter for a 
single flow and all packet lengths being equal to L. Then the total network throughput is 
given by 


e(8) — 


nBü - BL 
1+ D2, (860 - n2 (& +7.)) +0- Q — 8) = n8 = 72) T.) 
(7) 
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Since the denominator is bounded below by $77 , (sa Spy (& + By, it can be seen 
that 
TT E E 
OIX sa 


i.e., the total network throughput is bounded above by the reciprocal of the harmonic means 

of the physical bit rates of the n fiows. Thus, for example, if there are two flows with 

physical rates 2 Mbps and 4 Mbps then the total network throughput will be bounded by 
2 


TES 2.667 Mbps. Also, with equal packet lengths, we see that this total throughput is 


shared equally among all the flows. 


7 An Asymptotic Analysis 


We motivate a study of asymptotics with the plots in Figure 11. We notice that the fixed 
points appear to be converging, and there is not much variation in them for K > 15. In fact, 
for n < 20, there is little change in the fixed point even for K > 7. The reason for this is 
clear. The collision probability is small for these small values of n, and hence the retry limit 
K is rarely reached; increasing K beyond some value ceases to have any additional effect. 
Another inference we can make from these plots is that for large n, K = 7 gives a rather 
large collision probability. This may suggest that one should actually use a larger value of 
K in practice, say 10 or 15. This will not change the performance for small n and would 
improve the performance for large n. 

Thus we are motivated to analyse the fixed point for K — oo. Such an asymptotic 
analysis can be expected to yield closed form results that will provide information even for 
K as small as 15. A similar asymptotic analysis has also been carried out independently 
by Kwak et al in [7]; while their final results are the same as our Theorem 7.2, we have 
displayed an analytical form for the fixed point solution (see Theorem 7.1), and we derive 
our asymptotic results by taking a limit in this solution. Further, we also provide a relaxed 
fixed point iteration for computing the fixed point (see Section 7.1). 

To permit closed form analysis, let us take bọ = b slots, and by = p* x bo, where p > 1; 
hence, by Theorem 5.1, a unique fixed point still exists. The multiplicative increase is in 
any case a part of the IEEE 802.11 standard; we are generalising to an arbitrary multiplier 
in order to study the impact of the value of this multiplier. 

Assuming y < zi and taking K — oo, we see that 


l-py 
1-7 


GO) = x 


Note that the assumption that y < I does not affect the fixed point analysis presented 


earlier, since we will see in Theorem 7.2 that the fixed point in the limit K — oo is less than 
1 


Fh 
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Figure 11: The variation of the fixed points for K — 7 (top panel), K 
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= 10 (second panel), 


K = 15 (third panel), and K = 100 (bottom panel); b = 16, p = 2, and the number of nodes 
n = 2, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100. 
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Given y, G(y) is the probability of attempt of any node. Then using the batch Poisson 
version of the collision probability in Equation 3, the fixed point equation becomes 


yo ac dicm) " 


In order to obtain compact expressions, let us define 


|. n-—1 
n — bo 
A 
ze” 
x 
LambertW(x) 
—1 "p - 
z 
— (1/e) 


Figure 12: The LambertW function is the inverse function of ze”; notice that for x > -i, 
LambertW (x) < x, with equality only for x = 0. 


M P LambertW —1)e"?)—n5(p-1 
Theorem 7.1 The fied point is of the form y(n) = fan Ln 
Remark: For x > —+, LambertW (a) is defined as the inverse of the function ze”; see 
Figure 12. 


Proof: We proceed from Equation 8. Writing v = 1 — y, and using the definition of n, 


this equation can be rewritten as 


zd n(p—1) 
v=e Pe v 


Multiplying both sides by ņn(p — 1), we obtain 


n(p = 1) oS 


= 1 "nP — 
n(p — 1)e 7 


It follows from the definition of LambertW (-) that 
— 1 
mp) = LambertW (n(p — 1)e"?) 


from which the result follows by substituting 1 — y for v. M 
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7.1 A Relaxed Fixed Point Iteration 


The fixed point (7) can only be computed numerically. In this section we provide a relaxed 
fixed point iteration. With reference to Equation 8, let us define the function 


dietus a) 
and, with yo := A consider the sequence of values generated by the iterations 


Tei = (1 7 o) fk) + o (9) 


where 0 < a < 1. Notice that œ = 0 corresponds to the usual fixed point iteration, which 
will converge if f(y) is a contraction. The above iteration is called a relared fixed point 
iteration. We will now provide a condition on o that will ensure that the iterates converge 
to the fixed point. 

First of all, since f(y) is continuous, it is clear from the iteration in Equation 9 that if 
the sequence of iterates converge then they must converge to the fixed point. It is also clear 
that if, for each k, yx > f (yx) then the sequence {yx} is nonincreasing. This follows because 


Yeti = (1 — a) f yk) + ax S Ve 
if and only if 
fk) S Ye 


Thus, since x > 0, for the convergence of the sequence {7} it suffices to ensure that 


x > f(x) for all k. 
Now, it can easily be shown that the derivative of f(y) at y = I is given by D := 


n—l 


ale, and that, for yu 41 € Yk 


fna) — fex) € |D] (Ye — 0e) 


From this inequality we can see that to ensure f(x) < Yk, for all k, it is sufficient to ensure 
that, for all k, 


ID| (Yk — *&41) S Yer — fi) 


Using the iteration in Equation 9 this is equivalent to ensuring that, for all k, 


|D|(1 — o) — f(ye)) < aye — fn) 


Hence it suffices that 


IDI- a) <a 
or that 
ge Pl 
= [D|+1 
Thus, for example, with n — 10 nodes, bo — 16 slots, and p — 2, it relaxed fixed point 
iteration with a such that 222 < a < 1 will yield the unique fixed point y(n). 
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7.2 Taking n to oo 


We now wish to take n to oo and study the limit of the fixed point solution obtained in 
Theorem 7.1. For this we need the following properties of the LambertW function. The 
proof is provided in the Appendix. 


Lemma 7.1 
1. For a » 0, 


.. LambertW (axe) 
lim ——— — ———- 


1—00 x 


=1 (10) 


2. For 0 <a < 1, LambertW (aze?) € x 


3. For 0 < a < 1, the convergence in Equation 10 is from below. 


The following result is now obtained by applying Lemma 7.1 to the expression for y in 
Theorem 7.1. 


Theorem 7.2 


1. y(n) < i 


1 


3. limno nB 1 1n( 2—) 


p—1 


Remarks 7.1 


1. Theorem 7.2 provides explicit expressions for the collision probability and the fixed 
point for large K and a large number of nodes. We have already discussed above 
(recall Figure 11) that K only needs to be larger than 15 for the large K asymptotics 
to apply. Thus we see that for large n the collision probability is directly related to 
the back-off multiplier p, and is the reciprocal of this multiplier. 

2. We also see that ng, the mean attempt rate per slot, goes to In(ZEx), and hence the 
attempt probability per node (during back-off periods) behaves like O(+). This lends 
some support to the original assumption that from the point of view of a node the 
attempt process of the other nodes can be viewed as an independent process with i.i.d. 
batch Poisson arrivals in successive slots. 
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7.3 Asymptotic Aggregate Throughput 


Let us now consider n nodes handling n flows with all the flows having the same transmission 
rate, C. The aggregate throughput of the network is given by (compare with Equation 7) 


nde" L 


Q(B) = 1+ (apee (E F T.)) + (1 — 6775 npe P) Te) 


We infer from this equation that, as n — oo the aggregate throughput converges to 
0- iL 
1 lWL F 1 
seg + (0- 26 +T) + (a - 0-2) 1) 


The following result is then immediately obtained 


T(p) := 


Theorem 7.3 1. lim, 7(p) =0 
2. lim,_1 7(p) =0 
3. r(p) is maximised at 


Te 
Te+1 


p = ON a 
LambertW (-: - mu) + T 


€ 


Remarks 7.2 


1. The behaviour of the aggregate throughput as p goes to its two extremes is as expected. 
If p — 1 then the nodes do not increase their back-off intervals in response to collisions. 
The collision probability becomes large and the throughput drops to 0. Obviously, as 
p — oo collisions cause a drastic reduction in attempts essentially shutting the nodes 
off. 


2. In an attempt see what the above asymptotic results have to say about realistic network 
parameters, in Figure 13 we plot the aggregate throughput for finite K and finite m, 
using the formula in Equation 7 with equal transmission rate for all the flows. We see 
that the throughput increases steeply for 1 « p « 2, but is quite flat with p after p — 2. 
There is an optimal value of p, but unless p is very close to 1, the throughput is not 
very sensitive to p. It can be seen that the back-off multiplier used in the standard, 
i.e., p = 2, is adequate unless the number of nodes becomes very large. For Te = 17 
(slots), the third part of Theorem 7.3 returns p — 3.85, which compares well with the 
curve for n — 60 in Figure 13. 
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Figure 13: Aggregate throughput plotted vs. the back-off multiplier p for two values values 
of n. The network parameters are K = 10, bo = 16 slots, data packet length 8000 bits, 
packet overhead 592 bits, slot time 20 us, transmission rate for all flows 11 Mbps, fixed (rate 
independent) data packet transmission overhead 52 slots, collision overhead 17 slots. 


8 Comparison with a DTMC Model of the Back-Off Pro- 
cess 


The previous sections dealt with the following points: characterisation of the solution of 
the fixed point equation; analysis of system performance assuming that the system operates 
at the fixed point; asymptotic analysis of the fixed point as n — oo; numerical validation 
of the performance results obtained from the fixed point analysis in comparison with a 
discrete event simulation model. In this section we turn to the important issue of the 
relationship between an exact DT MC model of the coupled back-off processes and the fixed- 
point analysis. We proceed with the following assumptions: 


e The number of nodes n > 2 
e K is the retry limit 


e bo is the initial back-off, expressed as a multiple of the slot time. 


Exponential back-off with multiplier p > 1, i.e., bg = p'bo,1 € k € K. 


Back-off durations are geometrically distributed, or, equivalently, when a node is in 
back-off stage k, it attempts in the next slot with probability p 
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The n nodes all use the same back-off parameters. We focus only on the cumulative back-off 
time, since, as argued in Section 2.1, the evolution of the back-off process is what determines 
channel access to the nodes. With this in view, we index the slots during the system back-off 
time by t = 0,1,2,---. 

For the system with n nodes, at the beginning of slot t (more precisely at time t+), 
let, for 1 < j € n, s(t) (€ {0,1,2,---}) denote the back-off stage of node j; i.e., if 
s(t) = k, then node j will make an attempt at t + 1 with probability i If node j 
attempts at {+1 and there is a collision then st +1) = k+1, otherwise there is a success 
and Bia 1) = 0. If the node does not make an attempt at t+ 1 (with probability 1 — i) 


then st +1) = se) = k. Let S(t) denote the vector of the nodes’ back-off states 
at t. It is clear that S? (t), t = 0,1,2,---, is a Markov chain on (0,1, ---, K^. 

The dimensionality of the Markov chain S(?(t), t = 0,1,2,--- grows with n. It is 
therefore more convenient to work with the process that counts the number of nodes in each 
back-off stage. This will be a K +1 dimensional process for any number of nodes. To this 
end, given a state s € Z+” of the process S?) (t), define the number of nodes in the back-off 
stage k to be 


where J;.} denotes the indicator of the set {-}. Let M”) (t) denote the vector random process 


MPO, k € Z*. From the foregoing, it is clear that M(™ (t) is a Markov process taking 
values in the set 


K 
M) := {m : mx nonnegative integers; y Mk =n} 
k=0 


The following Lemmas will be used to establish that for bọ > 1, and p > 1, the DTMC 
M(”)(¢) is irreducible over the state space M). For two states m1 and mz in M? we 
write m; — m2 to denote that mə is reachable from mu, i.e., starting in state mə the 
DTMC reaches state m; in a finite number of steps with positive probability. 


Lemma 8.1 For bo > 1, and p > 1, for all m € MU), m — (n,0,0,---,0). 


Remark: This result asserts that from any state in M? the state in which all the nodes 
are in the initial back-off stage is reachable. 

Proof: Since bo > 1 and p > 1, for every k,1 € k € K, by > 1, i.e., the attempt 
probability of a node in any back-off stage is less than 1. Then starting from any state 
m € M?) the state (n,0,0,---,0) can be reached if exactly one node among the nodes in 
the stages 1,2,---, K attempts in m1 4- m2+---+mx successive slots, and hence succeeds. 
In each successive slot the number of nodes in back-off stage 0 increases by 1. This has 
positive probability. M 
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Lemma 8.2 For bo > 1, and p > 1, for all m € MU), (n,0,0,---,0) — m. 


Proof: See the Appendix. M 
The above two lemmas immediately yield the following result. 


Theorem 8.1 For bo > 1, and p > 1, the DTMC M?) (t) on M) is irreducible. Wl 


It follows that under the conditions bọ > 1 and p > 1, the DTMC M (t) is positive 
recurrent. Let x (^) denote the stationary probability measure on M™, 

For small values of K (e.g., 1 or 2) v?) can be numerically computed. Now given 2”), 
the collision probability y can be obtained as shown in Equation 11, where we have assumed 
that bj > 1 for all k. For a given m € M), in the next slot, the number of attempts by 
nodes in the kth back-off stage can be no more than mx. 


e (11) 
ES E (Qe€ed»eqec)gmne377)) 


) 
(meMG)j {am}, a; >2} x 
K 

) 


© 5 (se Yon (nt (02) G2 0-4) 
a;>1} : 


{meM(™)} tam, $0 
i= 


Denote such an attempt vector by a. Thus for given m, we have 0 < a € m. Then it 
is easily seen that in this expression the denominator is the total rate of attempts, and the 
numerator is the total rate of collisions. Note the terms x. ai, since we need to count all 
the attempts or collisions in each slot. 

À natural question to ask is whether the collision probability yielded by the above Markov 
chain is well approximated by the fixed point analysis. Sample results are shown in Table 1. 
Results are shown for K = 1 and K = 2. Two values of bọ are examined: bọ = 16, and also 
bg = 2. It can be seen that the fixed point analysis approximates the collision probability 
very well. 

A further exploration of the relationship between the exact Markov model and the fixed- 
point analysis is a subject of our ongoing research. 


9 Application to the Analysis of TCP Controlled File 


Transfers 
It is well known that over 9096 of the traffic in the Internet is generated by what are 
essentially file transfers between computers. All common applications such as e-mail, ftp, 


and web browsing basically involve no more than the transfer of files between one or more 
computers. Such file transfers are elastic in the sense that no guaranteed transfer rate 
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No. of DTMC FPA DTMC FPA 
Nodes | (K=1) | (K=1) | (K=2) | (K =2) 


Table 1: Collision probabilities obtained from the DTMC and from fixed point analysis 
(FPA) for K = 1 and K = 2. The table on the top is for bọ = 16 and the one at the bottom 
is for bo — 2. 
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is required, and, furthermore, the rate can vary during an ongoing transfer. 'Thus the 
transfer rates can adapt to the time-varying number of active connections in the network. 
In the Internet, such adaptation is carried out by the adaptive window mechanisms of TCP. 
When there is network congestion (indicated by packet loss, or specially marked bits in 
acknowledgement (ACK) packets) the TCP transmitter reduces its window, thus reducing 
the transfer rate. The default behaviour of the sender is to gradually increase the window, 
up to a maximum value, say, Wmax. In this section we are interested in using the analysis 
developed thus far in the paper to obtain formulas for the throughput of T'CP controlled file 
transfers rates over an IEEE 802.11 WLAN. 


9.1 Some Modeling Assumptions 


We will make the following assumptions: 


Al The files are infinitely long. Thus we do not deal with web transfers. Practically, this 
assumption means that our analysis applies to large file transfers, such as software, 
document, or media downloads. 


A2 The modulation scheme and bit rate of the physical connection between a pair of 
communicating wireless devices is ideally adapted (but fixed) so that there is no packet 
loss owing to bit errors. Further, the retransmission time-out at each TCP transmitter 
is large enough so that time-outs never takes place. 


A3 At the transmitter of each wireless device the capacity of the buffer is such that there 
is no packet loss. This assumption effectively holds in practice if the number of file 
transfer connections through a node is small enough so that the sum of the maximum 
TCP windows of all the connections is less than the buffer size. For, say, 10 connections, 
this would typically require a buffer of no more than 512 KB. 


A4 The file transfer throughputs are bottlenecked only by the rates they obtain over the 
WLAN. For example, the transfers could be between the wireless devices across an ad 
hoc WLAN, or, in the infrastructure case, between the wireless devices and devices 
attached to a high speed wired LAN to which the AP is attached. For transfers within 
a building or campus this assumption is practically valid since most wired LANs are 
based on 100 Mbps to 1 Gbps Ethernet. 


Owing to Assumption A1 it makes sense to talk about the long run time average through- 
put of a transfer. From Assumptions A2 and A3 it follows that the TCP window of each 
connection grows to its maximum value, and by Assumption A4, each data packet or ACK 
of all the TCP connections will be queued at the transmitter of one of the WLAN devices. 

Let us adopt the following connection model. There are m connections, indexed by 
j,1 <j € m. The source node of connection j is denoted by s(j), and the receiver node 
is denoted by r(j)(Z s(j)). Thus, for connection j, the TCP ACKs will queue up at the 
transmitter of node r(j). The data packet length for connection i is denoted by L; and the 
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ACK packet length by pe In general, each node will transmit data packets for some 
connections and ACK packets for other connections. 

In order to use the "saturated queues" analysis presented earlier in the paper, we make 
the following additional assumption 


A5 The configuration of the TCP connections and the sizes of their windows are such that 
the transmitter queues of the wireless devices never empty out. 


Remark: This assumption is made to permit us to use the fixed point analysis presented 
earlier in the paper. It, however, considerably restricts the scenarios to which the analysis 
wil apply. For example, the common situation of two or more devices simultaneously 
downloading files via an AP is not covered by our analysis. This is because the AP needs 
to send many more packets for each packet that each of the devices sends, and hence the 
device queues will empty out, violating our saturated queues assumption. 

We will utilise Assumption A5 as follows. Recall our discussion in Section 6. If all the 
n queues always have packets to send, then they always contend for the channel, and each 
successful attempt "belongs" to each of the queues with equal probability, Ł, 


9.2 A Formula for Connection Throughput 


Let us now focus only on the successful attempt instants. Such a success belongs to node 
i with probability l, The HOL packet at that node is then transmitted. If this packet 
is of length L and the transmission rate is C then a time Ł + T, elapses. If the packet 
transmitted is a data packet then possibly an ACK is inserted into the transmitter queue of 
the receiving node (note that if delayed ACKs are used then not every data packet causes an 
ACK to be generated). On the other hand, if the packet transmitted is an ACK packet then 
one or more packets are inserted into the transmitter queue of the receiving node. Thus the 
queues can be viewed as evolving only at successful polling instants. This is an important 
observation as it allows us to ignore the back-off periods while analysing the evolution of the 
packet queues. Note that this observation does not hold if there are finite rate open-loop 
arrival processes into the nodes, as these arrival processes will cause the queues to evolve 
even during back-off periods. 

From the above observations, we can now proceed by analysing the discrete time random 
polling model shown in Figure 14. The discrete “time” in this model evolves over packets. 
Note that we do not need to be concerned with packet lengths (data or ACK), or physical bit 
rates. We will see that all we need from this model is the fraction of polls to a queue that find 
packets of each type at the head-of-the-line. There are several TCP connections modeled as 
“chains” or classes of customers circulating between pairs of nodes. The populations of the 
chains are the TCP window sizes. If the delayed ACK threshold for a connection is greater 
than 1 (let us say 2), then at the receiving node for that connection 2 data packets give rise 
to one ACK packet. We can view this ACK packet as being a batch of 2, that is served 
together. 

The state of the random polling model is the position and type of each packet in each 
queue. This process evolves over packet times. It is easy to see that the evolution of 
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_--7TCP connection from node 2 to node 1 


random order 
of visit, with 
probability 1/n 


"es ~TCP connection from node n to node 2 


Figure 14: There are several TCP connections, modeled as "chains" of customers with a 
fixed population (the window size) circulating in a random polling network. The solid 
arrows between the queues show the direction of TCP data transfer for a connection, and 
the dashed arrows show the direction of TCP ACK transmission. The n transmitters are 
served in random order with equal probability for each node. 


this rather complicated process is Markovian. Analysis of this Markov chain will yield the 
following probabilities, that will be used in the throughput formulas. 


hij + the probability that at a polling instant the HOL packet at node à is a data packet 
from connection j 


REO : the probability that at a polling instant the HOL packet at node i is an ACK packet 
for connection j (for which node i is the receiver node, i.e., i = r(7)) 


By the observations made just before these definitions, we can conclude that the probabilities 
hi; and hi» do not depend on data and ACK packet lengths, nor on the physical bit rates of 
the connections. These probabilities will depend only on the maximum TCP window sizes, 
the delayed ACK thresholds, and the connection configurations (i.e., which nodes carry 
which connections). We also note that once we have these probabilities, the throughput of 
connection j can be immediately obtained as in Equation 12 (see also Equation 6), where 
(y, B) are obtained from the fixed-point analysis. 


UTE 


n—1 
BG — BN hs 


us (ack) 

X : n—1 X l Lj X : (able)? a n n—1 
1+ BG — B) hi j + hij —— | + To | +(-(1—-86) —»80-—95) ) Te 

= M) METEO) 

i- 


U:sG)-ij UwG)-i] 


(12) 


'This formula has the same form as the one in Equation 6. In the numerator the term 
B(1— 8)"— is the probability that node s(j) has a success, h,(;),; is the probability that the 
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HOL packet belongs to connection j, and when both these events occur connection j has 
a "reward" of L; bits. The denominator is the mean length of a back-off and and attempt 
cycle. 


Remarks 9.1 


1. To be technically correct Equation 12 should have been obtained as the ratio of two 
expectations with respect to the stationary distribution of the Markov chain describing 
the random polling model. We have shown only the final result in terms of the HOL 
probabilities at the polling instants, as this is simple and intuitively clear. 


2. In Equation 6 the HOL probabilities were obtained from the ratios of the open-loop 
arrival rates into the queues. In Equation 12, however, the HOL probabilities will need 
to be obtained from the packet level analysis of the random polling model shown in 
Figure 14. We will show how this is done in the next subsection. 


3. The denominator of the expression now includes a term for the service provided to 
TCP ACKs. 


4. We have used the fact that all data packets within TCP connection j have the same 
length L;, and the ACK packets within TCP connection j have the same size po 
If this were not the case then we would need to make a more elaborate definition of 
the HOL probabilities which would have to include the probability of finding packets 
of each possible length. 


9.3 Obtaining the HOL Probabilities 


Let A; be the throughput of connection j through its sender node s(j) in the random polling 
model shown in Figure 14. Thus A; the average number of packets of connection j that pass 
through the node s(j) per packet served in the polling model. 


Theorem 9.1 If at each success instant one of the nodes is polled with equal probability 
(i.e., we have the model in Figure 14) then he(j),j = Ajn. 


Proof: Let nsz); denote the fraction of packet services in the model of Figure 14 
during which the HOL position at node s(j) is occupied by a data packet of connection j. 
Since the mean time that a packet spends in the HOL position is n, by Little's Theorem we 
have 

Ta(j),5 = Aj 
Owing to random polling, the HOL position at node s(j) is observed by a Bernoulli process 
with probability of “success” equal to i, Hence by the result that Bernoulli “arrivals” see 
time averages, we can conclude that 


hs(5),5 = 15(5,5 = AWN 
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E 
Remarks 9.2 


1. If the throughput of ACKs for connection j through its receiver node r(j) is n 
then by the same argument as in Theorem 9.1 it follows that 


(ack) _ (ack) 
ITA Nn 


2. We note that the hypothesis of the Theorem 9.1 that “at each success instant one of 
the nodes is polled with equal probability” requires the saturation assumption, i.e., 
Assumption A5, to hold. There are TCP connection configurations for which this 
assumption will not hold. For example consider a single TCP connection from Node 1 
to Node 2. The TCP receiver uses a delayed ACK threshold of 2; i.e., it returns one 
ACK for two received data packets. Clearly over a large number of packet transmitted 
we cannot say that about half will come from Node 1 and the other half from Node 2. 
In this case the receiver node will tend to empty out and the saturation assumption 
will not apply. On the other hand if Node 2 was also sending to Node 1 then our 
analysis will apply. 


3. In view of Theorem 9.1 we need to analyse the random polling model and obtain the 
Ajs and the T and this will yield the HOL probabilities needed in the throughput 
formula. 


9.4 Comparison with ns2 Simulations 


In Figure 15 we compare the results from the analysis presented above and ns2 simulations; 
9596 confidence intervals are shown around the simulation results. For comments on the 
version of ns? used, see Section 5.2. The scenario simulated is that there are n nodes paired 
up with n other nodes; each node in the first group is performing a TCP controlled long file 
transfer to its corresponding member in the other group. The maximum receiver window 
for each TCP connection is 20 packets, the TCP packet length is 1 KB, and the receivers 
do not delay the ACKs (i.e., an ACK is returned for each received data packet). In this 
situation, of course, h;,; will be 1 whenever Node i is the source node of connection j, and 
heck will be 1 whenever Node i is the receiver node for connection j. The physical link 
rates are all 11 Mbps. The aggregate throughput over all the connections is plotted vs. the 
number of connections. We notice that the match between analysis and simulations is good, 
with the worst case error being about 6%. The simulation trace file showed that during the 
simulations there was no TCP time-out; thus our Assumption A2 held in this case. 
Another scenario that we evaluated was two nodes sending files to each other, simul- 
taneously. In this case the aggregate throughput predicted by the model is 2.4164 Mb/s, 
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Figure 15: Plot of aggregate TCP controlled file transfer throughput vs. number of simul- 
taneous transfers over an IEEE 802.11b network. Results obtained from the approximate 
analysis in the paper and from ns2 simulations are shown. We show 95% confidence inter- 
vals around the simulation results. In the ns2 simulation the IEEE 802.11 parameters are 
used; the data rate is 11 Mbps and the control rate (used to send RTS, CTS and ACK) is 
2 Mbps. 


while ns2 simulations return a 95% confidence interval of [2.4054,2.4153] Mb/s. In the sim- 
ulation, the throughput obtained by each transfer is approximately a half of the aggregate 
throughput. 


10 Summary 


In this paper we begin by providing a simplification and generalisation of a saturation 
throughput analysis of single cell IEEE 802.11 WLANs as provided by Bianchi [3]. We first 
focus on the back-off process which is analysed via a decoupling assumption and the solution 
of a fixed point equation. Our analysis provides a simple and general representation of the 
fixed point equation. The representation is insensitive to the distribution of the back-off 
times. We show that if the mean back-off durations for successive retrials are monotone 
nondecreasing then the fixed point equation has a unique solution. Then we provide general 
throughput formulas for open-loop arrival processes (e.g., UDP transfers). We recover the 
observation that connections with small physical rates dominate the throughputs of other 
connections. In fact, we show that in a simple case the network throughput is bounded by 
the reciprocal of the harmonic mean of the physical layer transmission rates. 
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We then turn to the special case of exponential back-off with an arbitrary positive mul- 
tiplier, p, and where we do not limit the number of retrials a node can make. This leads 
to simpler expressions which permit us to study the network performance as the number of 
nodes goes to infinity. For this case, we obtain a characterisation of the fixed point solution 
for the collision probability for each n. Then we take n to oo and obtain the limit of the 
collision probability and aggregate attempt rate that agree with the results of Kwak et al 
in [7]. We also provide a relaxed fixed point iteration for computing the fixed point for any 
finite n when the number of retrials is not limited. The asymptotic aggregate throughput is 
obtained and from this the optimal back-off multiplier p is also derived. 

For exponential back-off, and geometrically distributed back-off periods, the back-off 
process can be modeled via a discrete time Markov chain. We study this DTMC, and for 
some simple computable cases we compare the collision probability obtained from the DTMC 
with that obtained from the fixed point analysis. 

Finally, we show how the saturation throughput analysis can be used to obtain TCP 
controlled file transfer throughputs for some network scenarios. In this analysis we exploit 
the idea that for window controlled traffic the back-off process evolution can be decoupled 
from the packet service process, the latter being modeled by a random polling queue. 

Our analysis in this paper still leaves open several avenues for future work. Even for 
the saturated node case we have not taken care of the important physical phenomenon of 
capture which, owing to the use of direct sequence spread spectrum communication, could 
be significant for a small number of contending nodes. The fixed point analysis needs to be 
extended to nonhomogeneous nodes, i.e., the situation in which nodes use different back- 
off parameters, as would arise with IEEE 802.11e networks. Some extensions to the fixed 
point analysis in the homogeneous case and some results on the fixed point analysis for 
the nonhomogeneous case are provided in our companion paper [9]. Finally, the saturation 
assumption needs to be relaxed in order to permit the modeling of arbitrary open-loop 
arrival processes, or of general TCP transfer scenarios where the loads on the nodes are not 
balanced. 


11 Appendix 


Proof: (Lemma 5.1) 
We have 


l+y+7---+7* 


G (ee ee OT MUN RN 
) bo + ybi - 2b2 E + yh be +: + KK 


and we need to show that the derivative of this function with respect to y is negative. Taking 
the derivative we find that we need to show that 


K K K 
So bey D ED |< SC LS 5,360 
k=0 j=l k=0 


INRIA 


Fixed point analysis of IEEES02.11 35 


i.e., 
K K K K 
2,2 jh US Y gut m 
k=0 g=1 =0j=1 
or, equivalently, we need to show that 
2K min{n,K} 
Ee E j-w) 
n=1 j=max{(n—K),1} 
k-(n—3) 


min{n,K} 
amas (Ci KAFI 
k=(n—J) 


Now we consider each term 57. j(b; — bk) and show that it is non-negative. To 


this end, define 
m(n) = H0, k):j+k=n,1 > j,k < K}, 
where | - | denotes set cardinality. Then we can rewrite the previous expression as 


max{(n—K),1}+|9J-1 


((n — j) — 3) (ba; — bj) +n(bn — bo)lyien<k} 
j=max{(n—K),1} 


which is non-negative since, in the range of the sum, (n — j) — j > 0 and b,_; — b; > 0. M 
Proof: (Theorem 7.1) This is just a simple manipulation of the fixed point equation to 

get it into the form of the LambertW function. The fixed point equation is 
y=1- eoo) 


which can be rewritten as — 
]ip— 


(79) ere ons 


This expression can be rearranged as follows 


es | n(p—1) 
n(p — 1) e"? = TU 


7 


(1—) 
It follows, from the definition of the Lambert W function (and utilising the fact that p > 1) 
that (y) 
NWP — 
———— = LambertW e? 
— (n(p— 1) e") 


The result follows by rearranging the equation to extract y. Mi 
Proof: (Lemma 7.1) 


1. For x > 0, write z(z) = LambertW (aze*), i.e., 


z(x)e^? = aze” 
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It is easily seen that for x > 0, z(z) > 0, and z(z) f oo for z — oo. Now, taking 
natural logarithms, we obtain, for all x > 0, 


Inz(z) + z(z) 2lnaz + z 


Or 
z(x) naw +1 


x dnz(m) +1 


AL 


which, on taking x — oo yields the desired result since “2 and ao both go to 0. 


2. By definition, LambertW (xe*) = x, and LambertW is monotone increasing for positive 
arguments. Hence, for 0 < a € 1, LambertW (axe?) < x. 


3. Follows by combining the previous two parts. MI 


Proof: (Theorem 7.2) Observe that we can write LambertW (n(p — 1)e"?) as 


LambertW (2=npe"”), with pt being less than 1, by virtue of p > 1. The first two parts 
now follow upon using Lemma 7.1 in 


. LambertW (n(p — 1)e"?) — n(p — 1) 
V0) = — pambertWüp- Dew) 


The limit for n8 is obtained as follows. We have 
(l)e Cn 


Rearranging we have 


— In(1 — 
It follows that 
_ n 1 p 
mB = DB te pi) m 


Proof: (Lemma 8.2) Three cases arise depending on m 
(ii) There exists j', 1 € j' € K, such that mj, > 2 
(ii) There exists j',1 € j' € K, such that by m my —1 


Case (i): The sequence of transitions 
K K 
(n, 0,0, ---,0) est (mo, pnt mr, 0,---,0) "eek (mo, mi, X p- mi, 0,---, 0) > 
eoo (Mo, M1, Mo2,---,mx) has positive probability. The successive transitions correspond 
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to S jMk 2 2,1 € j € K, nodes attempting together and hence colliding, while the other 
nodes do not attempt. Each such transition results in the colliding nodes increasing their 
back-off stage by 1, until the desired state (mo, mi, m», --- , mx) is reached. 

Case (ii): This is the same as Case (i) except that we stop at the state 
(mo, M1, Mo, -- -, M 1,0, - 2 :,0). 

Case (iii): In this case m is of the form (mo, m4,m»,:--, mj —1,0,0, mı,0,0) where 
m, = 1 for some l, j' € | € K. There are two subcases here. 


(a) mo > 1 

(b) Mo = 0 

Subcase (iii-a): Consider the following transitions (n,0,0,---,0) — (mo — 1, > Mk + 

1,0,---,0) — (mo = 1, mı, Zo Me =F 1, 0,---,0) PE x (mo T: 1, M1, M2, >+, M1 DU 

2,0,---,0) — (mo — 1, mi, m»,--- , mj 4,0,0,---,0,2,0---,0) until the “2” reaches stage l. 

Then one of these nodes in stage l attempts yielding (mo, M1, ma,--- , mj —1, 0, 0, mu, 0,0). 
Subcase (ii-b): Define i = min(k : 1 < k € K, mx > 0] and consider the state 

(mi,migi, my —1, 0,0, mı, 0,0,0, ---,0), i.e., a left shift of m by i. Now consider the tran- 


sitions in Subcase (iii-a) until this state is reached. Then i more transitions in which all nodes 
attempt, and hence collide, will yield m = (0,0,---,0,0, mj, mi41,:-- mj ,0,0, mi, 0, 0). 
= 
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